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INTRODUCTION 


The  mathematical  formulation  of  the  most  general  problem 
of  Optimal  Control  can  be  considered  as  a  problem  of  Mayer 
subjected  to  unilateral  constraints,  i.e.,  to  certain  restric¬ 
tions  expressible  in  terms  of  inequalities  [ref.  2,  10,  12, 

17  and  19] . 

The  results  of  the  classical  calculus  of  variations  in 
their  usual  forms  cannot  give  a  general  solution  to  this 
problem  because,  among  other  things,  the  fundamental  relation 
of  the  calculus  of  variations,  i.e.,  the  equation  of  Euler- 
Lagrange,  is  Valid  only  in  the  case  of  points  interior  to 

the  set  of  admissible  points. 

To  date,  the  most  general  solution  to  this  problem  is 
given  by  the  Maximum  Principle  of  Pontryagin,  but  in  its 
present  form  this  principle  cannot  be  applied  in  certain 
situations,  and  its  validity  has  been  proved  in  particular 
cases  only  [ref.  2,  8,  15,  16  and  20].  It  is  our  intention 
to  give  a  derivation  of  this  principle  for  the  most  general 
case . 

This  derivation  corresponds  to  a  very  simple  and  very 
intuitive  geometrical  interpretation  in  the  event-  space, 
i.e. ,  the  state-space  with  an  extra  axis  for  the  time 
[ref.  22].  In  this  derivation  we  will  take  for  granted 
different  existence  theorems  when  the  geometrical  interpreta¬ 
tion  will  be  strong  enough  to  motivate  these  assumptions. 

The  existence  theorems  will  be  given  explicitly  in  another 

*  , 
paper  [ ref .  13  J  • 

With  the  help  of  this  method  we  will  also  treat  some 
examples  which  are  for  the  most  part  the  already  "classical" 
problems  of  Optimal  Control. 

*We  also  do  not  intend  to  give  in  this  report  a 
complete  analysis  of  the  connections  between  these  considera 
tions  and  other  parts  of  Mathematics,  as  Calculus  of  Varia¬ 
tions,  Theory  of  Semi-Groups,  etc.,  or  some  applications  to 
Theoretical  Physics,  Operations  Research,  etc.  All  these 
aspects  will  be  considered  later. 
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SECTION  I.  GENERAL  FORMULATION  OF  THE  PROBLEM 

We  assume  that  the  system  under  consideration  can  be 
completely  described  by  the  point  x  =  (x\  xn)  in  the 

n-dimensional  state-space  xn  and  that  its  evolution  is 
given  by  the  system  of  ordinary  differential  equations: 

x^  =  f1(x1,  .  .  .  ,  x11;  u1,  .  .  ur)  =  f 1  (x,u)  i=l,  •••>  n 

(1) 

The  vector  u  =  (u1,  . . . ,  ur )  in  Ur  is  called  the  control 
vector.  If  u  =  u(t)  is  known  and  if  appropriate  initial 
conditions  are  given  the  system  (l)  can  be  integrated  in  a 
unique  way  when  the  Lipschitz  conditions  are  satisfied. 

The  problem  under  consideration  is  to  select  a  partic¬ 
ular  vector  function  u  belonging  to  a  given  class  F  in 
order  to  meet  certain  requirements  which  will  be  discussed 
later  in  detail. 

The  specification  of  the  class  F  is  given  by  the 
particular  problem  involved.  This  specification  usually 
takes  the  following  form: 

u(t)  eF  if  and  only  if 

(i)  u(t)  is  piecewise  continuous 

(il)  gi(u)  >0  1  =  1,  . . . ,  k 

Usually  we  will  write  the  condition  (ii)  under  the  form 
ueft  .  For  instance  if  (ii)  is 


ft  will  be  the  unit  hypersphere  in  Ur,  or  if  (ii)  is 


:  1  , 


1  =  1, 


r 


-p 

ft  will  be  the  unit  hypercube  in  U  ,  etc. 

We  assume  that  the  initial  value  x(O)  =  and  the  final 

value  x(T)  =  of  the  solution  of  system  (l)  are  given. 

The  problem  is  to  find  a  vector  function  u(t)eP  such 

that : 

(i)  there  exists  a  T  >  0  for  which  the  integration  of 

(l)  with  the  control  u(t)  and  the  initial  condi¬ 
tion  =  x(0)  satisfies  x(T)  = 

(ii)  a  chosen  performance  criterion 


u(t )) 


dt 


0 


0 


(2) 


1) 


is  minimum. 

Some  particular  criteria  of  (ii)  will  be  discussed 
If  we  consider  the  case  f°(x,  u(t)  )  =  1  we  will  have 


i.e.,  we  will  require  the  process  to  take  place  in  the  minimum 
time . 


2)  In  the  more  general  case 

T 

fYx,  u(t  Mdt 

Jo  ^ 

can  for  instance  represent  the  total  energy  consumption, 
cost,  etc. 

3)  The  particular  problem  for  which  T  is  fixed  before¬ 
hand  can  be  treated  with  this  formulation  by  considering  a 

n* 

state  variable  x  for  which 


3 


1 


J 


and 


xn*(0)  =  0 


xn  (T)  =  T 


Let  us  return  to  the  most® general  case  (ii).  If  we 
Introduce  the  new  function 

f,t 

x°(t)  =  '  fYx,  u(t)^dt 

Jo  ^  ' 

we  can  add  the  differential  equation 

x°  =  fYx,  u(t)^ 

to  the  system  (l)  and  we  will  then  have  the  extended  system, 

x  =  f^x ,  u  ( t  ( 3  ) 

with 


X  = 

f^x. 

-P 

d 

X  = 

(x°. 

1 

X  y  •  •  •  y 

n  \ 

X  ) 

f  = 

(f°. 

r  y  •  •  •  y 

fn)  ' 

u  = 

(u1. 

p 

ro 

Va 

ur) 

We  assume  that  the  functions 


f^x,  u(t)^ 


i  =  0,  . . . ,  n 


are  defined  and  sufficiently  differentiable  for  all. 
(x,u)eXnUr 

Then  our  problem  can  be  stated: 

Find  a  function  u(t)eF  such  that 
(i)  there  exists  a  T  >  0  and  a  x°(T)  =  X°  for  which  the 
integration  of  equation  (3)  with  the  control  u(t)  and 
the  initial  condition  x(0)  =  (0,|1)  satisfies 
x(T)  =  (x° , 


pT 

[ ii )  X°  =  ffx,  u(tndt  is  minimum. 

virN  J 
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SECTION  II.  SET  OF  REACHABLE  EVENTS 


n+1 


Let  us  consider  the  (n+2)  dimensional  space  TX 

of  the  points  (t,  x°,  xt  •••"*>•  +1  deflned  by 

a  reachable  event  is  a  point  oi  1  x 

A  reacnaci  \  n  In  the  following  way 

a  function  u(t)eP  and  a  value  T  >  0  in  the 


t  =  T 


=  J  f^X,  u(t^ 


dt  i  =  0,  1 ,  • • • >  n 

subject  to  the  initial  conditions  *(0)  = 


The  set  of  all  reachable  points  will  be  called  R(?x) 
We  will  assume  this  set  to  be  dense  everywhere  and^its^ 
boundary  hypersurface  to  belong  ! 

this  boundary  hypersurface  S(^)- 


Example 


—  „  •  1  _  ,2  shall  move  with  the  condition 

A  system  x  =  u  snaxx  mwv 

1  2  2  2 

minimum  and  (u  )  +  (u  )  — 


'0 


this  means: 


r  •  o 

1 

j x 

=  u 

2 

u 

=  u 

u1)2 

p  2 

+  (u2)  <  1 

,0), 

i.e. ,  x°(0) 

and  x(0 )  =  (0,0),  i-e.,  x°(0)  -  0  and  x  (0)  0 


*  ^  r nnceot  in  Control  Theory  has 

The  application  of  this  co  PR  \  (private  communica- 
been  introduced  independently  by  E.  Roxin  IP*  x 
Lion)-  See  also  ref.  21- 
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Olic: 


The  set  of  reachable  events  )  is  nere 

points  (t,  x°,  x1)  such  that 


2  .2  o 

:°  +  x1  <  t2 


t  >  0 


in  other  words  R(0)  is  the  interior,  surface 
the  semicone  of  revolution  around  the  t  axis 
generatrix  x°  *  t .  The  boundary  hypersurface 
this  case  the  surface  of  this  semicone.  See 


included 
with  the 
S(0)  1 

Fig.  1. 


SECTION  III. 


PRINCIPLE  OF  OPTIMAL  EVOLUTION 


Let  us  call  0(6,, i2)  the  set  of  reachable  events 
satisfying  the  end  conditions.'  By  definition  <U6,,52) 


is  a  subset  of  RC^)* 
to 


More  formally  this  is  equivalent 


G(i 


’^2) 


P:PeR(|  )  and 


(x1 ,x2. 


n 


)p  -  6. 


} 


„  1-  ^<0  r  ( £  we  will  assume  that 

Among  the  elements  of 

A  is  the  one  with  the  smallest  x°.  We  will  assume  that 
such  a  point  exists,  in  other  words  that  the  infimum  of 
G(!  ,  £,  )  with  respect  to  x°  belongs  to 

1’By  definition  A  'is  the  end  point  of  the  optimal 
trajectory  and  the  function  u(t)  associated  with  this 
trajectory  is  the  solution  of  our  problem. 

We  see  immediately  that  by  definition  A  belongs  to 

the  boundary  hypersurface  S(^)  of  R^)- 


Theorem  I  „ ( t  \ 

Every  event  of  the  optimal  trajectory  belongs  to  S(y. 

The  proof,  by  contradiction,  is  immediate. 

This  theorem  constitutes  what  we  call  the  Principle  of 
Optimal  Evolution.  In  the  next  paragraph  we  will  express 
analytically  how  to  construct  trajectories- belonging  to 

S(^)- 

Examples 

!)  Let  us  first  consider  the  example  introduced  at  the  end 
of  Section  II.  In  that  case  G(0,a)  where  a  Is  a  given  value 
of  x1  will  be  the  set  of  all  points  to  the  right  of  the  rig 
branch  of  the  hyperbola  obtained  by  intersecting  the  semicone 
R(0)  by  the  plane  x1  =  a .  It  is  obvious  that  in  such  a  case 
there  is  no  point  of  G(0,a)  for  which  x°  is  minimum.  Then 
we  conclude  that  this  particular  problem  has  no  solution. 
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The  set  of  reachable  events  Is  defined  by 

(x°-t)2  +  (x1-  t)2  <  t2 
.  t  >  0 

l.e.j  a  semicone  ..entirely  situated  in  the  octan 

®  *" 

t  >  0  ,  x°  >  0  ,  x1  >  0  . 


G(0,l)  is  the  set  of  events  defined  by 
(x°-t)2  +  (1-t)2  <  t2 


t  >  0 


These  equations  describe  a  parabola  shown  in  Figure  3 • 

The  minimum  acceptable  value  for  x°  is  x  =0  and 
we  see  that  in  that  case 


T  =  1 

u1  (t )  =  -1 

u2 (t )  =  0 


Figure  3 


SECTION  IV.  GENERALIZED  HUYGENS '  PRINCIPLE 


In  geometrical  optics  there  is  a  simple  construction 
based  on  Huygens'  Principle,  which  gives  the  wavefront  at 
t  +  dt,  i.e.,  W(t  +  dt ) ,  when  the  wavefront  at  t,  i.e., 
W(t),  is  known  [ref. 11]. 


Figure  4 

From  every  point  of  W(t)  a  small  circle  of  radius 
cdt ,  called  wavelet.  Is  drawn  and  the  exterior  envelope  of 
the  small- circles  is  W(t  +  dt ) .  (See  figure  4). 

This  is  the  procedure  in  the  case  of  an  homogeneous 
isotropic  medium,  where  c  is  the  light  velocity.  This 


construction  can  be  generalized  when  these  wavelets  are  no 
more  circular,  for  instance  in  the  case  of  the  propagation  in 

a  homogeneous  anisotropic  medium. 

The  extension  to  the  case  of  wavelets  which  are  connected 

continuous,  and  differentiable  has  been  extensively  studied 
(Theory  of  Contact  Transformation,  Hamilton-Jacobi  Partial 
Differential  Equation,  etc.),  [Refs.  1,  4,  5,  6>  l4>  1  ,  22  . 

We  will  generalize  this  construction  to  the  case  of 
arbitrary  wavelets,  i.e.,  wavelets  for  which  the  conditions 
of  connectedness,  continuity  and  differentiability  have  been 

dropped . 


Examples 

1)  Let  us  assume  that  the  wavelet  associated  to  the  point 
is  the  line  segment  BC  (see  fig.  5)- 


Figure  5 

The  construction  of  the  new  wavefront  is  described  in  fig.  6. 
One  starts  with  a  wavefront  A  A'  A"  ....  and  one  obtains 
“the  new  wave frpnt  C  C'  C  ...  • 
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Figure  6 

2)  As  the  next  example  let  us  assume  that, the  wavelet 
associated  to  the  point  A  is  represented  by  two  points 
B  and  C  (see. fig.  7)- 


Figure  7 


The  construction  of  the  new  wavefront  is  described  in  fig.  8. 


w(t) 


Figure  8 

In  our  particular  problem  we  can  consider  the  inter¬ 
section  of  the  boundary  hypersurface  S(^)  by  an  hyperplane 
t  =  t  as  a  wavefront  W(t).  Then  the  wavefront  W(r  +  dt ) 
can  be  constructed  by  the  above  described  method.  In  fact 
this  will  allow  us  to  construct  the  whole  surface  S(^) 
from  the  point  (0,1-^)  which  is  the  intersection  of  S(^) 
by  the  hyperplane  t  =  0. 

Example 

In  the  second  example  of  Section  III,  W(x)  is  a  circle 
in  the  plane  (x^x1)  with  the  center  at  (t,t)  and  a  radius 
equal  to  t.  (See  fig.  9)- 


Figure  9 

The  wavelet  corresponding  to  a  point  A  will  be  a 
circle  of  radius  dt  and  of  center  (dt,dt)  relative  to 
(See  fig.  10). 


Figure  10 

The  construction  of  W(t  +  dt)  is  given  in  fig.  11 


t  t  +  dt 

Figure  11 


Usually  we  will  draw  the  wavelet  in  a  space  X  without 
the  factor  dt .  This  wavelet  is  then  the  mapping  of  ft  (set 
of  allowable  u)  into  Xn  by  the  relation 

x  =  f (x,u) 

The  example  under  discussion  is  described  by  the  differential 
equations 


u1  +  1 


u2  +  1 


with  uen  <=>  (u  )  +  (u  )  <  1 


1  2 

This  means  that  a  circle  in  the  u  u  plane  (fig.  12)  is 
mapped  into  a  circle  in  the  x°  x1  plane  (fig.  13)- 

In  general  the  wavelet  will  vary  with  x,  but  since  we 
assumed  that  the  functions  f(x,u)  are  sufficiently  differ¬ 
entiable  in  x  and  u,  we  can  nevertheless  be  sure  that  the 
successive  construction  of  the  wavefronts  W(t)  is  always 
possible . 
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SECTION  V.  ANALYTICAL  FORMULATION 
OF  THE  PRINCIPLE  OF  OPTIMAL  EVOLUTION 

Let  A  be  a  point  of  W(t)  where  W(t)  is  differentiable 
and  where  p  is  the  normal  to  W(t).  Then  the  question  arises, 
by  which  value  of  u  will  the  point  A  be  transferred  to  a 
point  B  of  W(t  +  dt)? 

If  we  define  H(x,p,u)  «  <p|f(x,u)>,  that  means  the  scalar 

product  of  p  and  f,  then  the  appropriate  control  function 

* 

u  is  determined  by 

u (x , p )  =  argmax  H(x,p,u)  y 

uen 

In  other  words  the  point  AeW(t)  will  be  transferred 
into  a  point  BeW(t  +  dt )  if  and  only  if  we  choose  the  control 
u  for  which  H(x,p,u)  is  maximum. 

Justification:  We  see  immediately  in  fig.  14  that 

Bew (t  +  dt)  is  the  point  corresponding  to  max  <p|f(x,u)>  and 
that  to  C,  not  on  W(t  +  dt )  corresponds  a  u*  such  that: 

<p|f(x,u*)>  <  <p|f(x,u)> 

For  the  time  being  we  will  assume  that  the  conditions  V 
determine  one  and  only  one  u. 

Theorem  II 

If  there  exists  a  control  u(t)  transferring  Aew(t) 
into  Bew(t  +  h)  and  CeW(t  +  h  +  k)  where  h  and  k  >  0, 
then  the  "topology"  of  W(t  +  h)  at  B  is  larger  or  equal 
to  the  "topology"  of  W(t)  at  A. 


We  define  the  symbol  "argmax"  by 

X  =  argmax  f(x)  if  and  only  if  f(X)  =  max  f(x) 

xe.B  x£s 
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V 


Figure  14 

Topology  of  a  wavefront  W  at  a  point  P  is  defined 
as  the  set  of  properties  of  W  in  the  neighborhood  of  P, 
for  instance:  existence,  continuity,  differentiability  of 
W  in  the  neighborhood  of  P. 

The  topology  at  B  is  said  to  be  larger  or  equal  to 
the  topology  at  A,  if  all  such  properties  of  W(t)  at  A 
are  also  properties  of  W(t  4-  h)  at  B. 

The  validity  of  this  theorem  is  easily  checked  in  the 
previous  examples.  The  general  proof  will  be  given  in  another 
paper.  The  practical  importance  of  this  theorem  will  be 
stressed  in  Section  XII  of  this  report. 

Corollary  II 

Let  W*(t  +  dt)  be  the  hypersurface  obtained  from 
W(t)  by  using  for  all  its  points  the  same  control  u  which 
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transfers  AeW(t)  into  BeW(t+  dt). 

By  definition  W*(t  +  dt )  and  W(t  +  dt)  intersect  at 
the  point  B,  and  by  application  of«  Theorem  II,  W*(t  +  dt ) 
and  W(t  +  dt )  are  even  tangent  at  the  point  B. 


In  particular  if  W(t)  admits  a  normal  tt(A)  at  A, 
so  will  ,W*(t  +  dt)  admit  a  normal  tt(b)  at  B,  by  defini¬ 
tion.  Moreover  the  Corollary  II  states  that  W(t  +  dt )  will 
also  have  a  normal  at  B  and  that  this  normal  will  also  be 


SECTION  VI.  GENERALIZED  HAMILTONIAN, 
FORMULATION  OF  THE  PROBLEM 


The  results  obtained,  in  the  previous  section  allow  us 
to  construct  a  trajectory  belonging  to  S(^)  if  we  know  the 
normal  p  to  W(t)  for  all  points  of  the  trajectory. 

For  this  reason  we  will  now  establish  a  system  of 
differential  equations  for  p. 

Let  us  define: 


<p  I  (  6x )  >  —  0  i  —  1  *  2 , 


<p  |  (6x)q>  =  -1 


>  VI 


where  the  (Sx^  are  n  independent  vectors  tangent  to 
W(t)  at  the  point  A  and  (6x)q  is  an  arbitrary  vector 
independent  of  the  (Sx^  and  directed  toward  the  inside  of 

W(t). 

We  will  assume  that  we  choose  the  same  control  vector 
for  the  point  x  and  for  all  the  n+1  points  x  +  (Sx)^ 
i  =  0,  1,  2,  . . . ,  n;  which  we  will  assume  to  be  in  the 
neighborhood  of  x. 

We  will  require  the  invariance,  along  a  trajectory,  of 

the  relations  VI,  in  virture  of  Corollary  II  of  the  preceding 

•  * 

section.  This  gives 


i  .  e  .  , 


<pj(6x)i>‘  =0  i  =  0,  1 ,  2 ,  .  .  .  ,  n 

<p  *  |  (  6x )  ±>  +  <P  |  (  6x )  1>  =0 


But 


(6x)^  =  A ( 6x) ± 


* 

time  . 


The  "dot"  indicates  differentiation  with  respect 


to 
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where 


CHAPTER  VII.  PONTRYAGIN'S  MAXIMUM  PRINCIPLE 

• 

The  results  obtained  at  the  end  of  the  previous  paragraph 
form  what  is  called  the  Maximum  Principle  of  Pontryagin. 

The  derivation  developed  in  this  report  gives  a  geo¬ 
metrical  interpretation  of  the  vector  function  p  intro¬ 
duced  in  the  formulation  of  this  principle  by  Pontryagin 
and  his  associates. 

The  interpretation  of  the  vector  function  p  as  the 
normal  to  the  boundary  hypersurface  S  (  £]_  )  of  the  set  of 
reachable  events  R(f|^)  allows  us  to  overcome,  with  the 
help  of  the  Theorem  II  previously  introduced,  most  of  the 
difficulties  which  usually  arise  in  this  application  of 
the  Maximum  Principle  to  a  particular  problem. 

These  are  generally  of  two  different  natures: 

1)  How  to  choose  the  initial  value  of  the  vector 
function  p,  i.e.,  p(0). 

2)  What  is  to  be  done  when  the  relation 

u  =  argmax  H(x,p,u) 
u££2 

does  not  determine  a  unique  value  for  u. 

With  our  geometrical  interpretation  it  is  easy  to  answer 
these  questions  as  we  will  show  in  the  section  devoted  to 
the  synthesis  of  the  solution. 
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SECTION  VIII.  CONTROL  WITH  NONZERO  INERTIA 


In  the  general  formulation  of  the  problem  we  have  defined 
the  class  P  of  acceptable  u(t)  by 

(i)  u(t)  is  piecewise  continuous. 

(ii)  g_,(u)_>0  1=1,  .  .  .  ,  k  i .  e  .  ,  uep 

For  some  particular  problems  the  restriction  (i)  is  not 
strong  enough:  the  rate  of  change  of  the  control  itself  is 
bounded,  for  instance  for  control  devices  with  nonzero 
inertia,  and  the  corresponding  reinforced  condition  (i)  should 
read: 

-X-  • 

(i)  u^(t)  is  continuous  and  u.  (t)  <  M 

Prom  the  theoretical  point  of  view  it  is  very  easy  to 
transform  a  problem  with  such  reinforced  conditions  (i)  into 
a  new  problem  with  usual  conditions  (i). 

Example 

If  tu(t)  is  to  be  continuous  with  u^t)  <  M  we  will 
replace  this  control  variable  u^(t)  by  a  new  one  called 
uk(t)  and  we  will  consider  ui(t)  as  a  new  state  variable 
to  which  will  be  associated  the  differential  equation 

ui(t)  =  uk(t) 
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SECTION  IX.  CHATTERING 


In  order  to  be  sure  that  the  solution  proposed  at  the 
end  of  the  paragraph  VI  satisfies  the  requirements  expressed 
in  the  general  formulation  of  the  problem  we  have  still  to 
check  that  the  control  function  obtained  by  this  method  be¬ 
longs  to  F. 

F  was  defined  as  the  set  of  functions  u(t)  such 

that 

(i)  u(t)  is  piecewise  continuous 

(ii)  uen 

The  condition  (ii)  is  satisfied  by  construction  but  it 
could  happen  that  the  condition  (i)  is  not. 

In  that  case  the  control  function  given  by  the  Maximum 
Principle  is  an  oscillating  function  of  infinite  frequency 
characterized  only  by  its  mean  value. 

When  such  a  case  occurs,  the  problem  has  no  solution, 
mathematically  speaking.  It  corresponds  to  a  variational 
problem  with  a  bounded  minimizing  sequence  whose  infimum 
does  not  belong  to  the  sequence. 

Example 

Find  the  shortest  continuous  and  differentiable  curve 
joining  two  given  points  A  and  B  of  a  plane  such  that  the 
tangents  to  the  curve  at  A  and  B  have  given  directions, 
one  of  which  at  least  being  different  from  AB  . 

From  a  physical  point  of  view,  nevertheless,  such  a 
pseudo  solution  is  not  without  interest:  it  gives  the  infimum 
of  x  (T),  i.e.,  a  value  impossible  to  actually  obtain  per¬ 
haps,  but  such  that  it  is  possible  to  get  closer  and  closer 
to  it  from  above'. 

In  a  specific  problem  when  such  a  case  arises  the 
problem  must  be  reformulated  and  the  modification  (i)*  given 
in  the  paragraph  VIII  is  to  be  performed. 
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SECTION  X.  MOVING  TARGET 


In  some  problems  the  terminal  condition  x(T)  =  is 

to  be  replaced  by  x(T)  =  (^(T)  where  ^(t)  is  a  given 
continuous  trajectory  and  T  is  the  time  when  the  target 
is  hit.  The  interpretation  in  the  space  of  events  of  such 
a  case  is  easy  and  we  see  immediately  that  the  Principle 
of  Optimal  Evolution  is  still  valid  with  all  its  consequences. 

In  such  a  case  T  is  the  smallest  value  of  t  for 
which  £  (T)  e  Since  we  assume  (^(t)  to  be  continuous 

we  can  even  say:  T  is  the  smallest  value  of  t  for  which' 
42(t)  e  S(^)  where  S  (f^ )  is  the  boundary  hypersurface 
of  R  (!,_). 


Example:  Brachistochrone  tracking  of  a  target  with  known 

trajectory. 

Let  S  be  the  system  to  be  controlled  and  T  the 

target.  We  assume  that  S  and  T  are  both  moving  in  the 

plane  x  x  .  We  know  the  trajectory  £(t)  of  the  target 

and  we  want  to  steer  the  system  in  order  to  hit  the  target 

in  minimum  time  (see  fig.  l6).  At  t  =  0  the  system  X  Is 

1  2 

assumed  to  be  at  the  point  x  =  0,  x  -  0.  The  speed  of  the 

o  2  •  2 

system  is  bounded:  x^  +  yfi  <  1.  Formally  the  problem  is: 

1  2 

Find  the  functions  u  (t)  and  u  (t)  which  give  the  minimum 


value 


for  which 


where 


x(T)  =  |(T) 

£(t)  is  a  given  function 

x(t)  is  the  solution  of  the  system 
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Figure  16 


•1 

1 

X  = 

u 

'2 

2 

x  = 

u 

with 

ul2 

+  u22  < 

and 

x(0) 

«  (0,0) 

1  2 

Solution:  let  us  consider  the  event  -  space  (t,  x  ,  x  )  where 

we  have  the  trajectory  ABCD  of  the  target  and  the  semicone 
of  the  reachable  events  of  the  system  (see  fig.  17). 


B  + 


Figure  17 


The  problem  will  have  a  solution  if  and  only  if  the 

trajectory  ABCD  intersects  the  semicone  of  reachable 

events.  Let  C  be  the  first  such  intersection,  we  see  that 

T  =  t  and  that  the  optimal  trajectory  for  the  system  is 
G 

OEC. 
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SECTION  XI. 


SYNTHESIS  OP  THE  SOLUTION 


Preliminary  Remarks 

The  Maximum  Principle  leads  to  a  solution  which  is  a 
relatively  strong  minimum  but  not  necessarily  an  absolute 
strong  minimum. 

Example : 

The  two  points  A  and  B  of  W(t)  (see  fig.  1  )>  not 
in  the  same  neighborhood  on  W(t),  can  be  transferred  to  the 
same  point  C  of  W(t  +  dt)  (see  fig.  19)- 


In  such  a  case  the  trajectory  t  passing  through  the 
point  D  for  instance  will  give  a  relative  strong  minimum 
but  not  an  absolute  strong  minimum. 


* 


Example 


with 


x1  =  u1 ( s in2  x1  +  sin2  x2 ) 

x2  =  u2 (sin2  x1  +  sin2  x2 ) 

,1.2  ,2.2 

uen  <=>  (u  )  +  (u  )  <1 
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In  the  following  developments  we  will  assume  that  such 
cases  will  not  arise . 

The  next  problem  is  to  analyze  the  situation  at  a  point 
A  of  W(t)  where  W(t)  is  not  differentiable. 


Figure  20 


to  which  corresponds  the  wavelet  of  fig.  21. 


A 


Figure  21 

In  such  a  case  the  maximum  principle  applied  to  p 
will  give  the  control  vector  corresponding  to  the  point  B, 
and  applied  to  will  give  the  control  vector  corresponding 

to  the  point  C. 

We  see  immediately  that  a  control  vector  corresponding 
to  any  point  D,  on  the  boundary  of  the  wavelet  between  B 
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and  C,  will  also  transfer  the  point  A  to  W(t+dt).  In 
such  a  case  we  see  that  the  transformation  of  W(t)  into 
W ( t  +  dt )  is  not  pointwise  one  to  one. 

Such  a  problem  could  also  arise  in  the  case  of  a  dif¬ 
ferentiable  wavefront  at  A  if  the  relation 

u  =  argmax  H(x,p,u) 
ueft 


does  not  determine  a  unique  control  u. 
Example:  Consider  the  wavefront  of  fig.  22 


md  wavelet  at  A  of  fig.  23> 


B 


Figure  23 


In  this  particular  case,  the  controls  corresponding 
respectively  to  B  and  C  will  both  transfer  A  into 
W (t  +  dt ) . 
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Conclusions 


If  we  consider  the  set  of  trajectories  obtained  by 
choosing  all  possible  values  of  p(0)  and  all  possible 
values  of  the  control  vector  satisfying  the  Maximum  Principle, 
we  will  obtain  a  generalized  field  of  trajectories,  corre¬ 
sponding  to  the  set  of  all  rays  in  geometrical  optics. 

This  generalized  field  may  contain 

1)  branching  points  (see  fig.  24) 


2)  disparition  points  (see  fig.  25) 


3)  indeference  regions 


i.e.,  a  subset  of  X  where  all  trajectories  passing 
through  all  its  points  and  corresponding  to  all  possible 
controls  belong  to  the  generalized  field. 

Example 
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SECTION  XII.  IMBEDDING  METHOD 


We  will  present  a  general  method  for  avoiding  the  dif¬ 
ficulties  arising  in  the  case  of  a  non-dif ferentiable  wave-  « 
front . 

In  the  formulation  given  earlier  the  wavefront  at  t  =  0, 
i.e.,  W(0)  is  reduced  to  a  single  point,  namely  (O,^). 

Let  us  now  assume  that  W(0)  is  a  n  dimensional  circular 
manifold  of  center  (0,^)  and  of  radius  e,  a  small  positive 
quantity o  The  orientation  of  this  W(0)  is  completely 
determined  by  its  normal  p(0).  By  applying  the  Maximum 
Principle  to  this  new  starting  wavefront  W(0)  we  may  obtain 
W(t)  for  any  t  >  0.  Because  of  the  Theorem  II  (see  Section 
V),  we  know  that  W(t)  is  differentiable  at  A,  a  point  of 
W(t)  corresponding  to  the  point  (0,^)  of  W(0). 

In  general  the  trajectory  (ies)  obtained  by  the  set  of 
points  of  W(t)  corresponding  to  (0,^)  of  W(0)  for  all 
t  >  0  will  belong  to  the  generalized  field  defined  in  the 
previous  paragraph  and  the  whole  field  will  be  the  set  of  all 
such  trajectories  cori’esponding  to  all  possible  choices  of 
W(0)  i.e.,  of  p(0).  This  method  gives  a  good  interpreta¬ 

tion  of  the  initial  choice  of  p(o)  and  its  application  will 
help  greatly  the  solution  of  most  practical  problems. 


*By  a  n  dimensional  circular  manifold  of  center  (0,^)  = 
(y^y1 ,  .  .  .  ,yri)  and  of  radius  e.  In  the  space  Xn+1  we 
mean  the  set  of  points  such  that 


n 


y  Pi(xi-y1)  =  0 


s 


i=0 

n 


,  I  1\.2  .  2 

(x  -y  )  <  e 


1=0 


where  the  vector  p  is  the  characteristic  orientation  of  the 
manifold.  Example:  if  n  =  2  the  circular  manifold  is  an 
ordinary  circular  disk  of  radius  e. 
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SECTION  XIII.  APPLICATIONS 


I.  A  second  order  system  subject  to  Inertial  forces  on_l£ 
shall  be  controlled  In  such  a  manner  that  a  fixed  target 
is  reached  In  minimum  time. 


This'  problem  Is  usually  given  in  the  form: 

How  to  choose  the  forcing  term  u  with  |u|  <  1 
to  find  the  quickest  path  of  integration  of 


in  order 


between  x 


dx 
dt2 


=  U 


dx 

and  x  =  a , 


=  b 


This  system  is  equivalent  to  a  system  described  by: 


with 


x1(0)  =  x2(0)  =  0 
x1(T)  =  a ,  x2(T)  =  b 


and  lul  <  1  <=>  uen  .  Q“ 

One  should  note  that  x°  =  1  and  x  (0)  =  0 
x°(t)  _  t.  It  allows  us  to  forget  the  variable  x 
analyze  the  problem  in  the  three  dimensional  space 


imply 

and 

t ,  x1 , 


Maximum_Principle 

H  =  <p\f>  =  P^1  +  P2f2  =  Pix2  +  P2U 

u  =  argmax  H  =  sgn  p ^ 
uen 


Differential  System_f or__P 


•  _  SH  _  n 

?1  Bx1 


dH 

p2  =  '  i?  *  Pl 


hence  p^  (t  )  = 

P2(t)  =  ~c1t  +  c2 

with  p(0)  =  (c1,c2)  i  °  e  .  ,  P1(0)  =  c1  and  Pg(0)  =  o£. 
Construction  of  the_trajectory_ correspond ing_to_a_particular_p(0);L 
We  have  found  u  =  sgn  p2 

Pg  =  -c^t  +  c2 

hence  , 

if  cx  =  0  and  c2  <  0,  we  will  have  P2(t)  <  0  for  t  >  0 

i.e.,  the  control  will  always  be  u  =  -1 

if  c  =  o  and  c2  >  0,  we  will  have  pg(t)  >  0  for  t  >  0 

i.e.,  the  control  will  always  be  u  =  +1 

C2 

if  c.  >  0  we  will  have  for  t  >  —  ,  Pp(t)  <  0 

1  °i  e- 

and  otherwise  ,  p  (t)  >  0 

c2 

i.e.,  the  control  will  be  u  =  -1  for  t  >  — - 

C1 

u  =  +1  otherwise . 

Cp 

if  c  '  <  0  we  will  have  for  t  >  +  —  ,  p  (t)  >  0 

1  C- 

and  otherwise  ,  p_(t)  <  0 

i.e.,  the  control  will  be  u  =  +1  for  t  >  +  — — 

.1 

u  =  -1  otherwise. 


% 
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In  other  words,  there  will  be  two  types  of  trajectories: 

-  one  type  will  be  formed  by  a  first  arc  with  u  =  +1 
followed  by  a  second  arc  with  u  =  -1 

-  the  second  type  will  be  formed  by  a  first  arc  with 
u  =  -1  followed  by  a  second  arc  with  u  =  +1 

with,  in  both  cases,  the  possibility  of  having  one  such 
arc  of  length  zero. 

The  integration  of  the  initial  system  for 
u  =  +1  '  gives  x"*‘(t)  =  Tj-  t2  +  x2(0)t  +  x^O) 

x2  (t )  =t  +  x2 (0) 

u  =  -1  gives  x1(t)  =  -  t2  +  x2(0)t  +  x^(0) 

x2  (t )  =  -t  +  x2 (0) 

Construction  of  the_wavefront _f or__t _=_T 

For  an  arbitrary  trajectory  of  the  first  type  we  will 

have 


+1 

for 

0  < 

t 

<  X 

with 

0  <  x  <  T 

-1 

for 

x  < 

t 

<  T 

This  gives  the  coordinates  of  a  point  reached  at  time  T 
after  starting  at  x2(0)  =  0  and  x2(0)  =  0: 

x1  (T )  =  \  t2  -  -|(T-x)2  +  t(T-t)  =  |  T2  -  (T-t) 


2^ 


x2(T)  =  t  -  (T-t)  =  2t-T 


For  a  given  T  the  system  a  represents  parametrically 
half  of  the  wavefront  W(T),  the  parameter  t  being  restricted 
to  0  <  t  <  T. 
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In  this  example  we  see  clearly  the  structureof  the 

1  2 

covering  of  the  whole  (x  x  )  plane  by  the  field  of  trajec¬ 
tories  from  the  origin.  Between  the  origin  and  an  arbitrary 
point  of  this  plane  there  is  one  and  only  one  such  trajectory, 
hence  the  solution  corresponding  to  the  given  end  point,  exists 
and  is  unique . 

If  we  consider  the  trajectory  OABC  we  can  see  how  the 
Theorem  II  applies:  At  any  point  X  on  OA  the  wavefront 
is  one  sided  differentiable  and  this  property  is  conserved 
for  any  point  of  the  trajectory  beyond  X.  At  any  point  Y 
on  AC  the  wavefront  admits  a  unique  normal  and  this  property 
is  conserved  for  any  point  of  the  trajectory  beyond  Y. 

A  direct  construction  of  the  field  of  trajectories  with¬ 
out  applying  the  Imbedding  Method  would  have  given  the  follow¬ 
ing  results: 

For  any  point  X  of  Otj/  we  have  the  choice  between 
u  =  +1  and  u  =  -1 ,  but  this  possibility  of  choice  is  con¬ 
served  if  and  only  if  we  choose  always  u  =  +1 . 

If  we  compare  this  result  with  the  content  of  Theorem  II 
we  can  say:  When  the  topology  on  a  wavefront  is  high,  the 
number  of  choices  is  low  and  reciprocal. 

II  Bushaw's  Problem,  i.e..  Problem  of  Minimum  Settling  Time 
for  a  Second  Order  System  with  Subcritical  Damping. 


*  o  . 
x  =1 


•1 

x 


-X 


X 


+x  +  u 


(1) 


For  simplicity  zero  damping  is  assumed. 


x1 (0 )  =  x2(0)  =  0 

with 

x1(T)  =  a,  x2(T)  =  b 


ko 


<3 

and  u  <  1  <=>  uefi  . 

One  should  note  that  x°  can  be  Ignored  as  In  I. 

Maximum  Principle 

- — —  — - — —  —  o 

h  =  <p|f>  =  Pjf1  +  p2f2  =  -plx2  +  p2^xl  +  u) 

Optimal  control:  u  =  argmax  H  =  sgn  p  (2) 

uen  d 


Figure  30 
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Analyt ical_ Forms  _ of _ the _Two_Types_ of _ Arcs 


The  integration  of  (l)  for  u  =  +1  gives 
x"*"(t)  =  rfcos(t+^)  +  1 

(5) 

'  x2(t)  =e^sin(t+^) 
with  eft  and  ip  determined  by  x(0). 

This  trajectory  is  a  circle  of  center  (+1,0)  and  of  frequency 
1  (c.p.s.)-  Similarly  the  integration  of  (l)  for  u  =  -1  _  . 
gives 

x1(t)  =  a/cos(t+f )  -  1 

(5)  ' 

x2  ( t )  =«^”cos(t  +p ) 

where  aft  and  ip  are  also  determined  by  x(0). 

This  trajectory  is  a  circle  of  center  (-1,0)  and  of  frequency 

1. 


Construction  of  the  Tra j ectory_Corresponding_to_a 
Particular  p(0) , _ i.e . , _to_a_Particular_0. 


We  have  found  u  =  sgn  p2 

..  and  P2  =  sin(t+0) 

hence 

u(t)  =  sgn  sin(t+0)  (6) 

In  other  words,  the  control  function  u  will  be  alternatively 
+1  and  -1  the  length  of  an  interval  being  7 r,  with  the  ex¬ 
ception.  of  the  first  and  the  last  intervals  which  may  have  a 
length  smaller  than  tr. 

Conclusions:  if  O<0<tt  u(t)  will  be 


+1 

on 

1 

•s 

O 

0] 

-1 

on 

[  7T—  0  , 

27T-0  ] 

(-l)n 

on 

[mr-0. 

(n+l)7r-0]  ,  etc.,  (see  fig.  31) 

+1  - 

u(t) 

1  JL 

0 

7T-0 

27T-0 

37 T-0 

4tt-0 

-1_ 

Figure  31 


Similarly  if 

-1 

+1 


-IT  <(})<  0 , 

on  [0 , 
on  [-0 


u (t )  will  be 

-0] 

-2+7 r] 


(-1) 


n 


An  example  for 


on  [-0+mr,  -0+(n+l)Tr] 

0  =  90°  is  given  in  fig . 
Switching  Curves 


32. 


If  we  call  N-l  the  locus  of  the  first  switching  of  a 
solution  starting  with  u  =  -1  (such  as  the  point  A  on 
fig-  32),  will  be  the  upper  half  unit  circle  around  (-1,0) 

(see  fig.  33)-  Similarly  P^,  the  locus  of  the  first  switch¬ 
ing  of  a  solution  starting  with  u  =  +1 ,  will  be  the  lower 
half  unit  circle  around  (+1,0). 
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Figure  33 


If  is  the  locus  of  the  i^*1  switching  of  a  solu¬ 

tion  starting  with  u  =  -1 ,  and  P  the  i1"*1  switching  of 
a  solution  starting  with  u  =  +1 ,  and  if  we  use  the  operator 
n  to  mean  an  inversion  through  (-1,0)  (see  dashed  lines  in 
fig.  32),  and  p  an  inversion  through  (+1,0),  we  obtain 
immediately: 


N 


2n 


pN 


2n-l 


N 


2n+l 


nN 


2n 


P2n  nP2n-l 


P2n+1  P?2n 


These  rules  and  the  knowledge  of  and  give  us  by 

induction  the  whole  switching  curve  [ref.  5]  (see  fig.  3^)  • 
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The  whole  set  of  trajectories  Is  shown  in  fig.  35- 


Figure  35 


Obviously  every  point  of  the  plane  is  reached  by  one 
and  only  by  one  trajectory. 
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Wave front  for  t  =  T 

Let  x(T,0)  be  the  value  of  x(T)  corresponding  to 
p(0)  =  (cos  0,  sin  0) 

(see  fig.  30). 

The  wavefront  for  t  =  T  Is  the  locus  of  x(T,0)  for  all 

0  >  1  •  e  •  y 

W(T)=<^X(T,0):  -7T  <_  0  <_  TTj> 

We  know  that 

u(t)  =  sgn  sln(t  4  0) 

By  Integration  of  the  Initial  system  with  such  a  control  we 
find  after  tedious  but  simple  computations  (see  Appendix) 
that : 


X(T,0)  =  x1 (T , 0 )  4  1x2(T,0)  = 


-2kei  (^+0)  _  e1T  4  el7rk 
If  0  <  0  <  tt(  1-k^ 

-2(k4l)  e1^)  -  e1T  -  e 


(7a) 


1 7rk 


If  7r(l-k1 )  <  0  < 


(7b) 


7T 


-2(k4l )  el(T+<^  +  e1T  4  elirk 
If 

1  (T+0 )  ,  „1T  _l7Tk 


-k17r  <  0  <  0 


-2ke^'* **  4  e1T  -  e 

If  -7T  <_  0  <_  -k.  7 r 


(7c) 

(7d) 


where  k-^  and  k  are  defined  by 


*  “ 

Please  note  that  1  =  V-l ,  the  Imaginary  unit,  and  not 
an  Integer  as  In  other  places.  » 

** 

By  definition  [x]  Is  the  largest  Integer  smaller  or 
equal  to  x. 
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This  general  structure  can  degenerate  and  give  the  follow¬ 
ing  situations: 

(l)  if  k  =  1  the  length  of  the  two  arcs  of  radius 
2 (k+1 )  is  zero  and  the  two  other  arcs  of  radius  2k  join 
tangentially  (in  fact  they  form  a  circle,  see  fig.  37) • 


Figure  37 

This  drawing  is  made  for  T  =  v  which  can  be  written 
T  =  O’tt  +  1*7 r  ,  i.e.,  k-^  =  0,  k  =  1  . 

(2)  if  k  =  0  the  two  arcs  of  radius  2k  disappear 
and  the  two  other  arcs  join,  but  not  tangentially  (see  fig.  38). 
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Figure  38 

This  drawing  is  made  for  T  =  ■?  which  can  be  written 
1  ^  1 

T  =  2  •  7T  +  0 ■  7F  ,  i.e.,  k1  =  2  and  k  =  0. 


Situation  (2)  happens  if  and  only  if  T  <  tt 

In  figure  39  we  have  drawn  the  successive  wavefronts 

for  T  =  tt,  -=2?-,  2tt,  -  *  •  and  the  trajectories 

3  3  s>  j>  ^  jj.  2 

corresponding  to  0  =  0,  +  ±  -g,  ±  ^ir ,  i r. 

If  we  consider  the  trajectory  OABCDEF  • • •  corresponding 

p 

to  0  =  —  —  7 r  we  can  repeat  the  remarks  of  example  I  on  the 
applications  of  Theorem  II  and  Corrolary  II.  In  short  at 
A  and  B  W(t)  Is  not  differentiable,  but  from  C  outward, 
W(t)  Is  and  remains  differentiable. 


APPENDIX 

DERIVATION  OF  THE  EQUATION  FOR  THE  WAVEFRONT  W(T) 

IN  BUSHAW'S  PROBLEM 

As  we  have  seen,  from  equation  5  in  Section  XIII,  the 
integration  of  the  system,  eq.  1  Section  XIII,  for  u  =  +1 
gives 

x1(t)=o$cos(t-t-'^)+l  (l) 

x2(t)  =a$sin(t+^) 

With  the  convention 

X  ( t  )  =  xX(t)  +  ix2(t)  (2) 

(1)  can  be  written 

X(t+tQ)  =  [X  (t  o )  -  l]eit  +  1  (3) 

■  and  in  the  particular  case  t  =  ir  we  have 

X  ( 7T+t  Q  )  =  -  X(to)  +  2  (4) 

Similarly  the  integration  of  the  system,  eq.(l) 
Section  XIII  for  u  =  -1  gives: 

x1(t)  =e^cos(t+^)  -  1 

o  ,  (5) 

x  (t)  =  a$; sin(t+V/) 

With  the  convention  (2),  (5)  can  be  written 

X(t+tQ)  =  [X(tQ)  +  1  ] eit  -  1  (6) 

and  in  the  particular  case  t  -  ir 

X  ( 7T+t  o )  =  -  X(to)  -  2  (7) 

The  formulae  (3)  and  (6)  can  be  replaced  by  the  unique 
formula 

X(t+tQ)  =  [X(tQ)  -  u]elt  +  u  (8) 

The  formula  (8)  is  valid  for  a  constant  u  which  can 
be  either  +1  or  -1. 

Similarly  (4)  and  (7)  can  be  replaced  by 
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(9) 


X(7T+to)  =  -  X(to)  +  2u 
We  have  seen.  Section  XIII,  eq.(6),  that 
u(t)  =  sgn  sin(t+0) 

In  Section  XIII,  eq.(8),  we  have  introduced 
decomposition 


with 

and 


T 

k 


77-(k] 


+  k) 


(10) 

the  unique 


In  other  words  k  is  integer 

and  0  <  kx  <  1 

In  the  following  computations  we  will  have  to  distinguish 
different  cases  depending  on  the  parity  of  k.  For  this 
purpose  we  will  introduce  further  decomposition 


with 


k  =  k2  +  2k^ 
k2=  0  or  1 


and 


This 


k3 

gives 


with 


integer. 

finally 

T  =  7 rfk-^  +  k2  +  2k^ ) 


0  <  k1  <  1 
k2=  0  or  1 
k^  integer 


(11) 


Let  X(T,$)  be  the  value  of  X(T)  corresponding  to 
p(0)  =  (cos  $ ,  sin  0). 

The  wavefront  for  t  =  T  is  the  locus  of  X(T,$)  for  all  0, 
i  .  e  .  , 


W(T) 


=  <^X(T,0):  -tt  <  0  <  TT^j  (12) 


53 


In  order  to  solve  our  problem  we  have  to  exhibit  a 


general 

expression 

for  X ( T , 0 )  for  all 

T  >  0 

and 

all 

-TT  <  0  < 

’  7 r  ■ 

We  will  distinguish  8  cases 

k2 

=  0 

and 

0  <  <p  <  tt( l-k1 ) 

Case 

no . 

1 

k2 

=  0 

and 

7r(l-k1)  <  0  <  TT 

Case 

no . 

2 

k2 

=  0 

and 

-k^7 T  <  0  _<  0 

Case 

no . 

3 

k2 

=  0 

and 

-TT  <  0  <  -k17T 

Case 

no . 

4 

k2 

=  1 

and 

0  <  0  <  7r(l-k1  ) 

Case 

no . 

5 

k2 

=  1 

and 

7r(l-k1)  <  0  <  TT 

Case 

no . 

6 

k2 

=  1 

and 

-k^7 T  <_  0  _<  0 

Case 

no . 

7 

k2 

=  1 

and 

-TT  <  0  <  -k17T 

Case 

no . 

8 

We 

see 

immediately  that 

X  (T  ,0  )  =  -  X(T,0+tt)  = 

-  X (T , 0- 

-it) 

(13) 

because 

W  (T )-  is 

symmetric  with  respect 

to  the 

center  of 

coordinates  X(0  ,<p)  =  0. 

With  the  help  of  relation  (13)  the  computation  of 
the  cases  3,4,7  and  8  will  be  particularly  easy  when  the 
results  of  the  cases  1,2,5  and  6  will  be  known. 

Case  no.  1 

=  0  and  0  _<  0  <  ir(l-k^) 

From  (10)  and  (ll)  we  see  that 


u  =  +1 

on 

o 

d 

i 

(I’l) 

u  =  -1 

on 

(  TT—  0 , 2  7T~  0  ) 

(1*2) 

1 — 1 
+ 

II  • 

d 

on 

(27T-0,37T-0) 

(1*3) 

u  =  -1 

on 

(  (2k^-l  ) tt~ 0 , 2 Trk ^ ~ 0 ) 

( 1 • 2k,  ) 

u  =  +1 

on 

(27rk2-0,T ) 

(l-2k3+l ) 
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Conclusion: 


in  the  case  no.  1  we  will  have  2k^+l  subarcs. 


By  assumption 


X(o,0)  =  0 


that  means  the  wave  propagation  starts  at  the  origin  of  the 
state-space.  By  application  of  (l°l)  and  (8) 

X(tt-0,0)  =  [X(O,0)-l]e1^Tr_0^  +  1 
-i 0  ,  , 

By  application  of  (1*2)  and  (9) 

X(2tt-0,0)  =  -  X(t r-0,0)  -  2 

-i  0  0 

=  -e  -  3 

By  application  of  (1*3)  and  (9) 

x(3 TT-0,0)  =  -  X(2tt-0,0)  +  2 

=  +e_i^  +  5 

etc . 

By  recurrent  applications  of  (l-i)  and  (9)  for  i  =  2,  3y-' 
2k^  we  obtain 


:(2-7rk3-0,0)  =  -  X(  (2k3~l  )tt-0,0)  -  2 


=  -e  -  4kq  +  1 


(proof  by  induction). 
By  application  of  (l-2k3+l)  and  (8) 


i  (T-27rk  +0 ) 


X(T,0)  =  [X(27rk3-0,0)  -  1  ]  e 

i(,jrk1+0)  irnk 
=  -  4k3e  1  -  e  1  - 

.  -ake1^)  -  e1T  +  eiTk 


(since  k„  =  0  in  this  case) 


Case  no.  2 


k2  =  0  and  7r(l-k^)  <  0  <  r 


Prom  (10)  and 

(ll)  we 

see  that 

u  =  +1 

on 

(O,TT-0) 

(2-1) 

u  =  -1 

on 

(  7T— 0  ,  0TT— 0  ) 

(2-2) 

u  =  +1 

on 

(27T-0,3Tr-0) 

(2-3) 

u  =  -1 

on 

(  (2^-1  )TT-0,2TTk^-0) 

(2-2k3) 

u  =  +1 

on 

(2k3~0,  (2k3+l)7T-0) 

( 2 • 2k3+l ) 

u  =  -1 

on 

((2k3+l)ir-0,T) 

(2-2k3+2) 

Hence  we  will 

have  2k ^  +2  subarcs  in  case  no. 

2. 

As  ( 2  ° i ) 

-  (1-i) 

for  1=1,  • • • ,  2k^ 

the 

computation 

in  case  no.  2 

will  be 

the  same  as  in  case  no. 

1 

up  to  the 

point 

X(27rk3~ 

■0,0)  =  -e10  -  4k3  +  1 

By  application 

of  (9) 

and  (2-2k3+l) 

X(  (2k^+l  )tt-0,0  )  =  -X(2Trk  3~0,0) 

+ 

2 

=  e  ^  +  4k  3  +  1 

By  application  of  (8)  and  (2‘2k3+2) 

i[T-(  (2k  +1  )tt-0)  ] 

X(T,0)  =  [X((2k,+1)7 r-0,0)  +  1  ]  e  J  -  1 

i(TTk  +0)  iTrk 

=  -2(2k^+l)e  1  -  e  1  -  1 

x  1  (T+0  )  IT  iTrk 
=  -2(k+l)e  -  e  -  e 
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Case  no.  3 


=  0  and  -k^7r  <  0  <  0 
We  use  the  relation 

X(T  ,0)  =  -X(T,0+tt)  =  -X(T  ,0*)  (13) 

X(T,0+7r)  is  given  by  case  no,  2  because  -k^  <  0  <  0 
implies  (l-k^)7r  <  0*  <  m  hence 

X(T,0)  =  -  [-2(k+l  )ei  (T+<£+70  _  e1T  -  ei7rk] 

-  +  eiT  +  ei1rk 


Case  no.  4 


k2  =  0  and  -tt  <  0  <  -k.  ir 

We  use  again  the  relation 

X(T,0)  =  -X(T,0+tt)  (13) 

X(T,0+tt)  is  given  by  case  no.  1  because  -ir  <  0  <  -k^Tr 
implies  0  <  0  +  m  <  (l-k^Tr,  hence 

X (T , 0 )  =  -  [-2kei  (T+<£+7r)  _  eiT  +  el7rk] 

=  -2kei(T+^  +  eiT  -  ei7rk 


Case  no.  5 


k2  =  1  and  0  <  0  <  Tr(l-k  ) 
Prom  (10 )  and  (ll)  we  see  that 
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u  = 

+1 

on 

(0,ir-(p) 

(5-1) 

u 

-1 

on 

(w-cp  ,2-ir-cp) 

(5-2) 

;  = 

+1 

on 

(27r-0,3Tr-$) 

(5-3) 

u  = 

-1 

on 

(  (2k  -1  )tt-0, 27^-0) 

(5°2k3) 

u  = 

+1 

on 

(27rk3-0,7r(2k3+l  )-0) 

(5*2k3+l) 

u  = 

-1 

on 

(ir(2k3+l)-0,T) 

(5-2k3+2) 

we 

will 

have 

2kg+2  subarcs  in  case  no. 

5.  As 

(5*i)  =  ( 2  - i )  for  i  =  1,  •••,  2k3+'l  fche  computation  in 
case  no.  5  will  be  the  same  as  in  case  no.  2  up  to  the  point 

X(  (2k^+l ) 7T—  0,0)  =  e  +  4k3+l 


By  application  of  (8) 

i[T-(  (2k  +1  )t r-0)  ] 

X(T,0)  =  [X((2k3+1)7 r-0,0)  +  1  ]e 

iirk„  i(Trk1+0) 

=  e  J  +  2(2k3+l)e  -  1  - 

,  -2ke1<T+«'>  -  e1T  +  e1* 


Case  no.  6 

k2  =  1  and  Tr(l-k1)  <  0  <  tt 


From  (10)  and 

(ID 

we  see  that 

u  =  +1 

on 

(0,7 r-0) 

(6-1) 

u  =  -1 

on 

(  If  —  0  y  271— 0  ) 

(6”2) 

1 — 1 

+ 

ll  • 

3 

on 

(2Tr-0,37T-0) 

(6-3) 

u  =  +1 

on 

(2Trk^-0,ir(2k3+l  )-0 ) 

(6’ 2k, +1 ) 

u  =  -1 

on 

(7r(2k,+l  )-0,ir(2k,+2)-0) 

( 6 • 2k,  +2 ) 

I — i 

+ 

II 

3 

on 

(ir(  2k^+2  )  -0  ,T  ) 

(6  •  2k,  +3) 

1 
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Hence  we  will  have  2k^+3  subarcs  in  case  no.  6.  As 
(6*i)  =  (2-i)  for  i  =  1,  2,  •••,  Sk^+l  the  computation  in 
case  no.  6  will  be  the  same  as  in  case  no.  2  up  to  the  point 

X(  (2k3+l ) tt-0 , 0 )  =  e-1<^  +  4kx  +  1 

By  application  of  (9) 


X((2k3+2)7r-0,0)  =  -X((2k  +1)7 r-0,0)  -  2 

=  -e~i<p  -  4k3  -  3 


By  application  of 
X(T ,0) 


(8) 

i 

=  [X(  (2k3+2)7r-0,0)  -  l]e 

i(k  +1)  i(Tr(k  +l)+0) 

=  -e  1  -  2(2k3+2)e  x  +  1 

-  -2(k+l)el(T+<>>  -  e1T  -  eilrk 


T- 


[7r(2k3+2)- 


0] 


+  1 


Case  no.  7 

kg  =  1  and  -k^Tr  <_  0  <  0 

X  (T , 0 )  =  -X(T,0+tt)  (13) 

but  X(T,0+7t)  is  given  by  case  no.  6  because  -k^ir  <  0  <  0 
implies  (l-k^)Tr  <  0+7 r  <  tt.  Hence, 

X(T,0)  =  -  [-2(k+l  )e1  (T+9+'+)_  eiT  _  giTrkj 

-  -2(lc+l)e1<T+*>  +  e1T  +  el7rk 


Case  no.  8 


kg  =  1  and 


-7 r  <  0  <  -k  7r 


X(T,0)  =  -X(T,0+tt) 


(13) 
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but  X  (T ,  0+7T )  is  given  by  ease  no.  5  because  -ir  <  0  <  -k  1?r 
implies  0  <  0  <_  (l-k^)TT.  Hence, 

X(T  ,0)  =  -  [ -2ke1  14+0+Tr]  _  glT  + 

=  +  e1T  -  eilrk 


These  results  can  be  summarized  by: 


Case 

no . 

1 

O' 

ii 

OJ 

ii 

and 

0  <  0  < 

7 T 

(1 

-kl} 

X(T  ,<fi)  = 

01  i(T40) 

-2ke  -  e 

IT 

4 

i7rk 
■  e 

Case 

no . 

2 

k2  =  0 

and 

7r(l-k^ 

)  : 

/< 

0  <  7T 

X(T,0)  = 

-2(k+l ) 

ei(T40) 

- 

e 

IT 

-  e 

i7rk 

Case 

no . 

__3 

o 

II 

OJ 

and 

_k]_7r  < 

0 

< 

0 

X(T,0)  = 

-2(k+l ) 

ei(T40) 

4 

e 

iT  , 

4  e 

i7rk 

Case 

no . 

_4 

o 

II 

OJ 

ii 

and 

VI 

b 

i 

< 

- 

k^rn 

X(T,0)  = 

01  i  ( T -f  0  ) 

-2ke  v  ^  H-  e 

IT 

- 

i7Tk 

e 

Case 

no . 

i — 1 

II 

OJ 

ii 

and 

o  <  0  : 

<  7 r( 

l-k1) 

X(T,0)  = 

01  i(T40) 

-2ke  v  -  e 

IT 

4 

l7Tk 

e 

Case 

no . 

_6 

kg  =  1 

and 

ir(l-k 

)  < 

0  <  7 r 

X(T,0)  = 

-2(k4l) 

gi  ( T+0 ) 

- 

e 

IT 

-  e 

i7rk 

Case 

no. 

_ V 

k2  =  1 

and 

-k-^TT  < 

0 

< 

0 

X(T,0)  = 

-2 (k4l ) 

gl  (T40  ) 

4 

e 

iT  i7rk 

,  e 

Case 

no . 

_8 

kg  =  1 

and 

-7T  <  0 

< 

- 

k-^7T 

X(T+0)  = 

01  i(T40)  , 

-2ke  v  4  e 

IT 

_ 

i7rk 

e 
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We  see  immediately  that  the  results  in  cases  7>  >  6,  7,  and 
8  are  identical  with  the  results  in  cases  1,  2,  3 ,  and  4 
respectively .  In  other  words,  has  no  influence  on  the 

structure  of  the  final  answers.  Therefore,  we  will  thus 
restrict  our  attention  to  the  cases  1,  2,  3,  and  4  only, 
in  the  consideration  of  the  shape  of  the  wavefront. 


Case  no.  1 


0  <  0  <  tt^I-^ ) 


v  m  \  i(T+0)  iT  ,  imk 

X(T+$)  =  -2ke  v  ^ -  e  +  e 


This  is  an  arc  of  a  circle  of  radius  2k  and  of  center 

/  IT  ,  imk.. 

( -e  +  e  .  ) . 


Extremal  values: 


X(T,0)  =  -(2k+l)e1T  +  ei7rk 


(l4) 


Tangent : 


X(T,7r(l-k1))  =  (2k+l)el7rk  -  e1T  (15) 


X(T,0)  =  -2kie 


i(T +0) 


Extremal  tangents:  X(T,0)  =  -2kie' 


(16) 


X(T,7r(l-k1) )  =  +2kie’ 


(17) 


)ase  no.  2 


7r(l-k1  )  <  0  <  7r 


X(T  ,<p)  =  -2(k+l)e1  -  e1T  -  el7rk 


This  is  an  arc  of  a  circle  of  radius  2(k+l)  and  of  center 

/  IT  imk-, 

(-e  -  e  ). 


Extremal  values: 


X(T,7r(l-k1))  =  (2k+l)el7rk  -  e1T  (18) 
X(T,vr)  =  (2k+l)eiT  -  e17rk  (19) 


The  dot  on  X(T ,0)  means  a  derivation  with  respect  to  0 


X(T,0)  =  ^  x (t ,0 ) 
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Tangent : 


X(T,0)  =  -2 (k+1 ) ie 


i(T+0) 


Extremal  tangents:  X(T,7r(l-k,  ))  =  2(k+l)le 


Imk 


X(T,tt)  =  2  (k+1 )  ie 


IT 


(20) 

(21) 


Case  no.  3 


-k^m  <  0  <  0 

X(T,0)  =  -2(k+l  )e1  ('T+^/l  4  elT  4  e1^ 


This  Is  an  arc  or  a  circle  of  radius  2(k4l)  and  of  center 

/  iT  ,  i7Tks 
(e  4  e  )  . 


"1  1<C  TT  "1  T* 

Extremal  values:  X(T,-k^ir)  =  -(2k4l)e  4  e 


X(T,  0)  =  -(2k4l)elT  4  e17rk 


Tangent : 


X(T  ,cp)  =  -2(k4l)ie 


i  (T40 ) 


(22) 

(23) 


Extremal  tangents:  X(T,-k^7r)  =  -2(k4l)ie 


ikm 


X(T,0)  =  -2 (k4l ) ie 


IT 


(24) 

(25) 


Case  no.  4 


-TT  <_  0  <  — k-^TT 


X  ( T ,  0  )  =  -2ke1  (T+t4  )  +  e1T 


iirk 


This  is  an  arc  of  a  circle  of  radius  2k  and  of  center 
,  IT  iirk>. 

(+e  -  e  )  . 


Extremal  values: 


Tangent : 


X(T,-tt)  =  4(2k4l)elT  -  e17rk 
X(T,-k1  ir)  =  -(2k4l  )elirk  4  e1T 
X(T,0)  =  -2kie1(T+!^ 


(26) 

(27) 


Extremal  tangents:  X(T,-tt)  =  2kie 


IT 


X  (T  ,  -k^ir )  =  -2kie17Tk 


(28) 

(29) 
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By  comparing  the  formulae  (l4)  to  (29)  we  can  check 
the  continuity  of  X(T ,<p)  and  the  continuity,  up  to  a  real 
factor,  of  X(T,0),  i.e.,the  continuity  of  the  direction  of 
the  tangent. 
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